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The electronic structure and magnetic properties of the Jahn- Teller-distorted 
perovskite KAgF3 have been investigated using the full-potential linerized aug- 
mented plane-wave method. It is found that KAgF3 exhibits significant quasi-one- 
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dicular to the z axis) and J (along the z axis) about 0.04, although the sublattice of 

magnetic ion is three-dimensional. The strong quasi-one-dimensional antiferromag- 
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I. INTRODUCTION 

Low- dimensional magnets have unique electronic and magnetic properties. For example, 
the physical properties of quasi-two-dimensional (2D) high-temperature superconductor-^, 
spin-Peierls compound 3 -, Haldane chain 4 and spin-ladder 5 as well as spin-frustrated 
systems^ have been extensively investigated. In most of the low-dimensional magnets, the 
magnetic cation sublattices consist of planes or chains which are kept reasonably far apart. 
Materials which exhibit three-dimensional (3D) magnetic ion sublattice and low-dimensional 
magnetic behavior have been of interest to both experimentalist and theorists.- - — Up to date, 
the electronic orbital ordering (00) is considered as an essential factor in determining the 
strong spacial exchange anisotropy in these materials. 

KAgF 3 is a material which exhibit strong spacial exchange anisotropy within 3D magnetic 
Ag 2+ ion sublattice. 13 It was first synthesized in 1971 by Odenthal et alM and received, 
as well as other silver fluorides, attention due to the pursuit of superconductivity^— in 
transition-metal compounds other than the cuprates. Previous experimental studies^ 3 -*^ 
have provided the evidence of its antiferromagnetism with TV=64 K as well as the insulating 
characteristic below this temperature. In contrast to the undoped cuprates that are quasi-2D 
antiferromagnetic (AFM) insulator, KAgF 3 was reported as a quasi-one-dimensional (ID) 
antiferromagnet using LSDA+ U method.— However, more efforts are still needed to well 
understand the origin of the quasi-lD antiferromagnetism. 

In the present work, we investigate the electronic structure and the magnetic properties 
of KAgF 3 using density functional theory (DFT). The calculated J±/\ J|=0.04 confirms that 
the compound is a quasi-lD antiferromagnet. Here, J (J±) refers to exchange constant 
along (perpendicular to) the z axis. The quasi- ID antiferromagnetism can be understood 
by the C-antiferro-distortive 00 which can be obtained in the presence of on-site Coulomb 
repulsion U. The picture of 00 and AFM insulating state is thus very similar to that of an 
isoelectronic and isostructural compound, KCuF 3 . 

II. CALCULATION METHODS 

Our electronic structure calculations are performed within the full-potential linearized 
augmented plane- wave framework. 1 - The generalized gradient approximation (GGA) of 



Perdew-Burke-Ernzerhof form is adopted.' 19 To include the on site Coulomb interaction 
GGA+U (t/=on-site Coulomb repulsion strength) approach is used with U e ff—U-J (J 
is the exchange interaction) instead of U. 2Q Here in our work, on-site Coulomb repulsion 
U is applied to Ag 4d orbitals only. The Muffin-tin sphere radii are chosen to be 2.24, 
2.09, and f.85 bohr for K, Ag, and F atoms, respectively. Within the Muffin-tin sphere 
the electrons behave as they were in the free atom, and the wave functions are expanded 
using radial functions (solutions to the radial part of Schroinger equation) times spherical 
harmonics. Out of the Muffin-tin spheres wave functions are expanded using plane waves. 
The value of R.MT^max (the smallest muffin-tin radius multiplied by the maximum k value 
in the expansion of plane waves in the basis set) is set to 7.0. We use 500 k points (i.e., 
192 k points in the irreducible wedge of the Brillouin zone) for the integration over the Bril- 
louin zone. Self-consistency was considered to be achieved when the total energy difference 
between succeeding iterations is less than 10~ 5 Ry/unit cell. The present setup ensures a 
sufficient accuracy of the calculations. 

In our study, we carry out total energy calculations with three distinct spin states of 
(1) ferromagnetic (FM), (2) yl-AFM (ferromagnetism in xy plane, AFM stacking), and (3) 
G-AFM (antiferromagnetsm in xy plane, AFM stacking). 

III. RESULTS AND DISCUSSIONS 

The calculations were performed for orthorhombic structure proposed by Mazej et al. 13 , 
space group Prima (No. 62), with the lattice constants a=6.2689 A, 6=8.3015 A, c=6.1844 A 
(see Fig. 1). In our calculations we use for each AgF 6 octahedron such corresponding local 
xyz coordinates that the z axis is along the b axis and the x and y axes are approximately 
along the inplane nearest neighbor Ag-Ag directions. In KAgF3, the Ag ion sublattice forms 
a pseudocubic structure, in which Ag-Ag distance along z direction is slightly shorter than 
that in the xy plane. There are two types of F sites. One connects the Ag ions along 
the z-axis (named as F ap ) and the other connects the Ag ions in the xy plane (named as 
F p i). The AgF 4 planes are slightly puckered since the F p / ions have a 0.31 A displacement 
from the plane formed by Ag ions. Each Ag and its six near F ions form a Jahn- Teller 
(JT) distorted AgFg octahedron. The distortion leads to the alternating long and short 
Ag-F bonds along the x and y axes. The four shorter Ag-F bonds are almost identical. The 



cooperative JT distortion provides a signature of an orbital ordering (00) at the Ag 2+ sites, 
which determines the electronic structure as well as the magnetic properties. 

In order to clarify the electronic structure and the origin of the quasi- ID magnetism 
in KAgFs, our calculations are designed in two stages. We start the study within GGA. 
To account for a possible lattice distortion, we carried out an optimization of the atomic 
positions, keeping the unit-cell parameters fixed and relaxing the atomic coordinates. The 
structural optimization shows that the Ag-F ap bond is elongated by 0.012 A and the in-plane 
longer (shorter) Ag-F pi bond is shortened (elongated) by 0.085 A (0.038 A) as compared 
with the experimentally observed structure. As a result, the large in-plane distortion of 
about 0.3 A found experimentally, decreases to 0.2 A theoretically. Moreover, the GGA 
results for the optimized structure show that the FM, ^4-AFM and C-AFM spin states are 
not stable and converge to nonmagnetic (NM) metallic solution. Fig. 2 shows the total 
and the orbital-resolved density of states (DOS) for the NM KAgF 3 . The distortive AgF 6 
octahedron crystal field splits the Ag 2+ Ad 9 orbitals into fully occupied ti g (xz, xy, yz) 
and three fourth occupied e g (z 2 -y 2 , 3x 2 -r 2 ) orbitals. Both the anti-bonding Ag d z 2_ y 2 and 
ds x 2_ r 2 bands cross the Fermi level with large band widths of about 3 eV and have strong 
covalency^ with the F orbitals. The band center of the d z 2_ x 2 is higher in energy than that 
of d 3y 2_ r 2 orbital. Thus, the JT distortion generates a crystal- field splitting which is not 
large enough to separate the two e g bands. Even the GGA calculations for the experimental 
structure still give a NM metallic solution for all the enforced spin states, in accordance 
with the previously published work.— Unlike the electronic structure of GGA optimized 
structure with two type of e g bands crossing the Fermi level, the distortion alone separates 
the e g bands with only one of them cross the Fermi level (not shown). Therefore, the GGA 
calculations can not reproduce the insulating AFM nature of KAgF 3 . This is in contrast to 
another fuoroargentate Cs2AgF 4 in which JT distortion alone can lead to the orbital ordered 
insulating state.— - — 

As the second stage, we concern the correlation interaction by GGA+ U approach. In 
general, strong on-site Coulomb repulsion U increases the localization of d electrons and 
favors orbital ordering and/or magnetic order. However, the magnitude of the correlation 
for the Ag 2+ ion is not known. In order to determine its value, we performed a structural 
optimization as a function of U (1-6 eV) and found the U value with which the equilibrium 
structure matches the experimental one. This scheme to extract the Coulomb parameter is 



feasible because the on-site Coulomb repulsion U determines the orbital polarization of the 
e g states for a large extent, which, in turn, causes the structural JT lattice distortion. The 
optimized inequivalent Ag-F distances d\, d,2 and d% are shown in Fig. 3. We find that the 
optimized structure for [7=5.0 eV is the closest to the one obtained experimentally. And 
thus this U is adopted as the most reliable one. This value is smaller than the commonly 
used f/=7-9 eV value for Cu 2+ but still reasonable, because the 4d orbitals of an Ag 2+ ion 
are less contracted than the 3d orbitals of a Cu 2+ ion. 

The GGA+ U (5 eV) total energy calculations of various spin states shown in table I 
indicate that the ground spin state is the ,4-AFM insulating state. The obtained ground 
spin state is in accordance with the previous LSDA+U results. 13 This indicates that the 
on-site Coulomb repulsion is important in determining the AFM insulating properties in 
KAgF 3 . We note that the difference in total energy between the ^4-AFM and G-AFM spin 
states is as small as 23 meV/4f.u., while the energy difference between the ^4-AFM and FM 
spin states is about 272 meV/4f.u.. We should notice the fact that ^4-AFM and G-AFM spin 
states have opposite ordering in the xy plane but the same ordering along the z direction, 
whereas both the ^4-AFM and FM spin states have the same ordering in the xy plane but 
different ordering along the z axis. This reveals that AFM exchange interactions along the 
z direction are more favorable and robust than FM coupling in the xy plane, indicating a 
strong spacial exchange anisotropy in KAgF 3 . 

To gain additional insight into the spacial exchange anisotropy, we evaluate the mag- 
netic exchange constants along and perpendicular to the z axis numerically in terms of the 
Heisenberg spin Hamiltonian: 

ff = J^S i .S i + J L J]S fc -S« (1) 

i,j k,l 

where the first (second) term refers to summing over all nearest neighbors along (perpen- 
dicular to) the z axis. By mapping the obtained total energies for each magnetic state to 
the next nearest Heisenberg model, the exchange interactions J and J± are 

J = 8$2( Efm ~ E A-afm) ( 2 ) 

J± = T&p( E A-AFM - E G-AFm) ( 3 ) 

With the spin 5=1/2 of a 4d 9 Ag 2+ ion, we get J=136 meV and Jj_=-5.8 meV, indicating 
strong AFM coupling along the z axis and a much weaker FM coupling in the xy plane. 
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|J^|/J=0.04 reflects strong quasi-lD magnetism in KAgF 3 . Thus, it has been concluded 
that KAgF 3 is another material exhibiting quasi-lD magnetic behavior within 3D magnetic 
ions sublattice, which is in accordance with the results of Mazej et al..— The relatively 
large J and J± values are consistent with the findings that the DFT electronic structure 
calculations generally overestimate the magnitude of spin exchange interactions.— 

The quasi- ID magnetic behavior can be understood from the ground state electronic 
structure and 00 of KAgF 3 . The total and the orbital-resolved DOS for the v4-AFM spin 
state experimental structure are shown in Fig. 4. The Coulomb repulsion U pushes the 
occupied and unoccupied 4d levels downward and upward respectively and opens up an 
insulating gap. It is obvious that the one hole states mainly occupy the higher level d z 2_ y 2 
bands and the orbital polarization of e g states is enhanced. Actually, due to the strong 
Ag-F covalency^ the one hole spreads over the six fluorine atoms of the AgF 6 octahedron. 
This is also reflected by the magnitude of local spin magnetic moments within each muffin- 
tin sphere, ±0.62 //e/Ag, 0/xs/F ap and ±0.09/ib/F^ (see table I). The hole states have 
alternating d z 2_ y 2 and d z 2_ x 2 symmetry for Ag 2+ ions in the xy plane (see the last panel in 
Fig. 5) because of the cooperative JT distortion. The in-plane d z 2_ x 2J ' d z 2_ y 2 00 repeats 
along the z-axis and C-antiferro-distortive 00 is formed, which is clearly shown in Fig. 
5. The overlap of Ag-F-Ag along the z-axis (see the upper and middle panel in Fig. 5) is 
larger than that in the xy plane (the lower panel in Fig. 5) leading to the dominate magnetic 
interaction along z-axis. The C-antiferro-distortive 00 would immediately give ^4-AFM spin 
state, according to Goodenough-Kanamori- Anderson rules. Also, our GGA+ U calculation 
results confirm this mechanism since it shows that the ^4-AFM state is indeed more stable 
than the C-AFM and FM states (see Table I). Therefore, the results of our calculations 
have clearly illustrated that KAgF 3 is an orbitally ordered quasi-lD AFM insulator, in close 
analogy to KCuF 3 . 

IV. CONCLUSIONS 

By GGA and GGA+ U electronic structure calculations, we have found the quasi-one- 
dimensional (ID) antiferromagnetism of KAgF 3 with J±/\ J|=0.04, which may stimulate 
further experimental studies. Here, J (J±) refers to exchange constant along (perpendicular 
to) the z axis. The quasi- ID antiferromagnetism can be understood by the C-antiferro- 
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distortive orbital ordering. The orbital ordered insulating state can only be obtained when 
the on-site Coulomb repulsion of Ag 4d electrons is in consideration. With U=5 eV for Ag 
4d electrons the optimized structure is in the best agreement with experiment. The picture 
of orbital ordering and AFM insulating state is thus very similar to that of an isoelectronic 
and isostructural compound, KCUF3. 
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FIG. 1: (color online) Crystal structure of KAgF3. The blue, gray and yellow spheres represent K, 
Ag and F atoms, respectively. 

FIG. 2: (color online). The total and partial DOS of GGA optimized structure for nonmagnetic 
state obtained in GGA. The Fermi level is set at zero energy. 

FIG. 3: (color online). Deviation of relaxed inequivalent Ag-F distances d±, di and d% from the 
experiment as a function of the Coulomb repulsion U. The data of U=0 denotes the GGA optimized 
results. 

FIG. 4: (color online). The total and partial DOS of KAgF 3 for ,4-AFM spin state from GGA+U 
with U=5 eV. The up and down panels in each DOS plot denote the spin-up and spin-down states 
respectively. 

FIG. 5: (color online). The contour plot of spin density in xz(upper panel), yz (middle panel) and 
xy (lower panel) planes of KAgF 3 for yl-AFM state from GGA+U (5 eV). The solid (dot) lines 
depict the spin-up ( spin-down) states. 

TABLE I: Electronic structure of KAgF 3 in FM, ,4-AFM and G-AFM spin states obtained by 
GGA+U with U=5 eV. The total energy difference (AE, meV/4f.u.), the energy gap (E g , eV), 
the local spin moment (SM, fig) of each Ag, apical F (F ap ), and planar F (F p j) atoms are shown. 
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23 
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±0.6 





±0.09 



10 



Fig. 1 Zhang.eps 
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Fig. 2 Zhang.eps 
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Fig. 3 Zhang.eps 
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Fig. 4 Zhang.eps 
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Fig. 5 Zhang.eps 
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